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Recently, Rindler and Ishak have argued that the bending of light is, in principle, changed by the 

presence of a cosmological constant since one must consider not only the null geodesic equation, but 

also the process of measurement. I agree with the fact that both must be considered. Here, on the 

basis of the mathematically exact solution to the classical bending problem, and independent of the 

cosmological constant, I clarify the approximate argument found in the vast majority of texts (new 

and old) for the measured value of the bending of light for a single source and show that the result is 

f~^ ' in part due to the almost perfect cancelation of two terms, one of which is seldom considered. When 

("^ , one considers two sources, this cancelation is of no consequence, and, for example, if the sources are 

f^ ■ opposite with the same associated apsidal distance, the approximate argument gives the rigorously 

CNJ ' correct answer (up to numerical evaluation), an answer which is unaffected by the presence of a 

►^ , cosmological constant. 

)_^ ■ PACS numbers: 04.20.-q, 04.20.Cv, 04.20.Ha 

(N; 

O ! I. INTRODUCTION 

trk' Recently, Rindler and Ishak [l| have corrected a long-standing error in the literature concerning the cosmological 

I i' constant (A) and the bending of light for a single source, an error perpetuated by the author [^. Here I reexamine 

the classic subject of the bending of light (with a concentration on a solar mass deflector - not gravitational lensing 

\N ' in general) and arrive at some further results which should be of wide - spread interest. The background geometry is 

J^ r the spherical vacuum given by 

^1 ds^ = ^ + r^de^ + sin^ edcj,^) - f{r)dt\ (1) 

where 



and we restrict our analysis to 9 = -n 12. 



o : /(o = 1 - ^ — ^ (2) 



X 



II. NULL GEODESICS 

The non - radial null geodesic equation associated with ([T]) given ([2]) can be written in the form 

(S)' -<- K - "' + 2"' (3) 

dcp 

where u = ni/r and the maximal value of u is assumed fixed and given by Uj, = m/r^ < 1/3 irrespective of A [3|. 
Equation ^ is solved exactly, up to sign, by [1| 




'/'W = J u^.W •^|2\/^^,^,\/i^l, (4) 

where 

e^{l-2u^){l + 6u^)>l, (5) 



so that 



and 



The function l{u) is given by [5| 



A = 6u^-1 + VQ>6u^> 0, 

B = 6u^-1-Ve < 2{3u^ - 1) < 0, 

k{u) = 4w + 2uj, - 1 

Ve + k{u) > 2u^ > 0, 

Ve-k{u) > 2(1 -3uj,) > 0. 

l{u) = -2ul + Wj, - 2u^u + li - 2w^ > 



(6) 

(7) 
(8) 

(9) 

(10) 

(11) 



and T is the incomplete elliptic integral of the first kind of purely imaginary modulus Jg . The solution ^ is shown 
in Figure [H Note that there is an apse at B where u = u^ and <t>{u^) = V Uj, < 1/3 Q- 



For u — Q we have 



*<°) = ^>/^^ VtV^ 



Along the orbit, Up is distinguished by the condition 



TT _ / 9- k'^jup) I l u^ -Up A 

2 "V -i{up)B -^ rV A 'Vs 



(12) 



(13) 



Table |T] gives a numerical summary for grazing incidence with the Sun based on (accurate) numerical approximation 
to the exact solution ^. 




FIG. 1: Part of the null geodesic A^ given by (|4]) in the u — (f) subspace of ^ with ([2]). The vertical axis is exaggerated for 
clarity. The line OB is the line of apsides so the complete trajectory N is obtained by adding the reflection of the part of A'^ 
shown about OB. A is explained in Table U Exhibition of the entire orbit, without use of asymptotics, derives from the use 
of u instead of r. 



It is fair to say that the vast majority of texts (new and old) consider only equation ([3]), in the equivalent form 

cPu 



u = 3u , 



(14) 



by way of a first order approximation. This procedure gives rise to the classic total deflection 4w_ (for A = 0) which, 

as Table H] shows, is an excellent approximation to the total deflection defined by 2A. (Indeed, the second order 
approximation, given by Au^+u'^{'^^ — 4), is very much better still! See also Appendix A.) However, we do not live 
at M = 0, or even at Up (in the inner Oort cloud) but at m_ which has to be considered far away from w = in the sense 
that 0(0) — (/)(u„) ~ 959". 62. Presentations that rely on considerations of ([3]) (or (fT4| ) alone are, therefore, misleading 
[W^ . Yet, "the bending of light" ~ 1".75 has been measured [11|, and without going to m = ! To understand this it 
is necessary to give an operational definition to "the bending of light" which we now consider. 



Quantity (units) 

0(0) (rad) 

20(0) - TT EE 2A (") 

4«o (") 

0(^t©) (rad) 
20(^3,) - ^ (") 



TABLE I: u^ = M, 



iable ^. 



Value ~ 

1.57080 rad 

1".75051 

2A-0". 722465 lO"'^ 

2A-0". 336958 10"^° 

p = 1096.38 

1.56615 rad 

-1917".49 



III. THE BENDING OF LIGHT 



The "bending of light" {ip) is defined here to be the combination of the deflection prior to the apse plus the deflection 
after the apse, 



Prior to the apse we have 



■0 = A„ + Aout- 



Ain = (piUe) 



(15) 



(16) 



where e stands for the emitter and 4>{ue) follows from Q. If lie = 0, and for grazing incidence with the Sun, Ai„ = A. 
Approximations to Ai„ are carried out in Appendix A. 

The deflection Ao„t, unlike Ai„, involves the process of observation and the difference between two angles: Xj the 
angle between the tangent to the null geodesic A^ and the direction of the deflector as measured by a timelike observer 
on an orbit of constant u in the u — 4> subspace of ([T]) (the measured position of a star relative to the center of the Sun 
during a total solar eclipse in the classic bending of light experiment), and 0, the angle deflned between the position 
of the deflector and the perpendicular to the line of apsides. This relates the position of the star ^ 6 months after the 
eclipse as explained below. The angles are shown in Figure [5] along with their evolution along N for various observers. 
The relationship of (j) to theory and observation is explained in Figure [S] 




FIG. 2: Demonstration of the angles x and 4>. The evolution of these angles for various observers along A'^ is also shown. The 
angle x shrinks monotonically along A'' from the apse, where x = ''"/S, to O where x = 0. The angle shrinks monotonically 
along N from the apse, where (j> — it/2, to p where = and then increases to A at O. This is discussed further below. 



eD 
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FIG. 3: The undeflected position of a star is defined to be tliat position determined by radial null geodesies from the star. In 
the Figure AB is the line of apsides. The perpendicular becomes a radial null geodesic on rotating the Sun through n — (j) about 
the Earth, or equivalently, the Earth about the Sun by tt/2 + 0(itg,). We have assumed no parallax. The theoretical value for 

(j> is then 7r/2 — (j}{ua.)- In practice, in the classical experiment, the undeflected positions are used to form a template and the 
time for alignment is ~ 182.4 days (let us say 6 months). 



We have 



A{u)out = x{u) ~ 4>{u) = x(u) + • 



w-5. 



(17) 



where (j){u) is given by (g]). It is clear from Figure[2]that A{u^ )out — 0, A{up)out — x(up) and A{0)out — </'(0)— 7r/2 = A 
(for a solar deflector). We are interested in A{u^)out- 

It is important to note that the reference point chosen for the definition of the measured angle x does not in fact 
change A{u)out- For some other reference point we have 



A{u)out = x*(w) + a - 4){u) = x*iu) + a + cj){u) 



IT 

2' 



(18) 



where the angles are shown in Figure ID The angle x*(u) is measured, 4'{u) follows from Q and a from the geometry 
as shown. In all cases X = X* + Q^- These angles are measured and calculated in the 9 — 7r/2, t = const subspace of 
([T]) with ^ in which we use the standard formula 



cos(t) 



d^S' 



for the angle t between the two directions defined by d* and (5*. 

From ([!]) with ([2|) and p9ll we find the theoretical prediction for x [li 

sin(x) = 



(19) 




(20) 



where we have written x for xi^)^ f for /(") ^^'^ fs f*-"^ /(^s)- (Because we have chosen a timelike observer on an 
orbit of constant u to measure Xj we have the restriction u > u-^. This does not change in any significant way the 
current presentation.) Approximations to x are carried out in Appendix B. Now if, as shown in Figure SI we were to 
choose x*j a similar calculation gives 



cos X ) = — ' 

U \ ./f /lJi-y2 4- f _ 1 



^^(^r + f- 



(21) 




FIG. 4: Deflection in the r — plane showing the variation of the reference point chosen for the definition of the measured 
angle x- As an example, for x the reference point is the center of the sun but for x* the reference point is the edge of the solar 
disc. Both Ai„ (shown in insert) and Aout simply measure the deviation from the horizontal (straight - line motion). 



with 



cos(a) = 






(^)^+/-l 



(22) 



Form ([^ . (PT|) . and ([^ it can be shown that X — X* + ct- Now we could choose the reference direction to be the 
measured position of the deflected star in which case x* — ^^'^ Qfj in numerical value, = x- With the foregoing 
understood, we emphasize the choice x over x* just for simplicity in what follows. First, however, it is instructive to 
look at the weak - field limit / ^ 1. In this limit wc obtain x = arcsin(rE/r), x* = ^nd a = arcsin(rE/r) as we 
would expect. 



From ([U, P?]) and (PO)) . with all currently acceptable values of A, and for grazing incidence with the Sun, wc find 



m 



A(wg^)<,„,:^l".75049 (-). 



This result is best understood in the following way: 

-<^(ue) = 0( 
and 






-1917".49 (-) 



-1919". 24 (-). 



(23) 



(24) 



(25) 



It is, of course, crucial here that we evaluate (piu^) from ^. The almost exact cancelation by addition of these two 
contributions at w„, along with the assumption that u,, — 0, gives the famous result '>jj{u^) ~ 1".75. This cancelation 
is, of course, not unique to the Earth [ly|. Table HIl demonstrates this. 

The vast majority of texts (new and old) consider only the deflection of a single source and argue that the observed 
bending of light derives from the fact that we can set both the emitter and observer at u = and so 



i) = 2A„j = 2(t){u ~ 0) - vr = 2A ~ l".75. 



(26) 



Whereas this argument arrives at a good answer, on the basis of the exact solution to ([3]) I have shown that in fact 
holds and so these approximate arguments leading to (|26p do not adequately explain the measured deflection of 



TABLE 11: M^ = 



u^, u variable [14 



Distance (AU) 

1/214.75 

1/214 

1/213 

1/212 

1/210 

1/100 

1/10 

Mercury [15] 

Venus 

Earth 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

50 

500 

1096.38 

5000 



-630397 

-609289 

-592452 

-579625 

-559334 

-199621 

-19197.5 

-4956.37 

-2651.59 

-1917.49 

-1257.87 

-367.030 

-198.545 

-98.0566 

-62.0036 

-36.6341 

-2.08795 



1.36666 



2X 

630398 

609290 

592452 

579626 

559335 

199623 

19199.2 

4958.12 

2653.34 

1919.24 

1259.62 

368.781 

200.295 

99.8071 

63.7541 

38.3846 

3.83846 

1.75051 

0.383846 



2Ao.ut 

0.07467 

0.16402 

0.23500 

0.28881 

0.37335 

1.54953 

1.74862 

1.75038 

1.75048 

1.75049 

1.75050 

1.75051 

1.75051 

1.75051 

1.75051 

1.75051 

1.75051 

1.75051 

1.75051 



light at the Earth when a single source is observed. If we set Ue = 0, so that Ai„ = A, we must still account for 
AoMt and for that I have shown that the measured value is the result of the almost perfect cancelation of two terms, 
p4)) and ((25l) . one of which is seldom considered, namely ([25]) . It is also this term that formally introduces A by way 
of (|20|) and ^. This cancelation is discussed in Appendix C. If one chooses not to use the center of the sun as the 
reference point for the measurement of Xi this cancelation remains. It is then the cancelation oi x* + a with that 
explains measured deflection of light. 



IV. WIDE-ANGLE DEFLECTIONS 



Nowadays the classical bending of light test is extended over the entire sky and is used with interferometric methods 
at radio wavelengths to put (remarkable) limits on the PPN parameter 7 ^T\ . This is done on the basis of the Shapiro 
- Ward formula Il8|, which in our notation reads as 



*~2u^(H-cos(xH)). 



(27) 



Note, however, that in (|27p x{'^) is often taken as an observed quantity, and not derived from (|20p . Here we compare 
this relation with ^j: as defined in p^ . The situation is summarized in Table Hill where larger values of the apsidal 
distance are considered up to x == 7r/2 and the observer is at u = u„. Whereas there is no measurable difference in 
the two definitions (-0 and '^), they are not equivalent and the cancelation discussed in the previous section plays the 
dominant role. This is also discussed in Appendix C. 



V. DIFFERENTIAL MEASUREMENTS 

In practice it is the measurement of differential deflections in opposition that is also of importance |19| . Consider 
two sources, 1 and 2, which have associated with them the same u^ (but not necessarily for grazing incidence with the 
Sun) that are seen on opposite sides of the deflector (in the classic case, during a total eclipse). Then the differential 
deflection is given by 



(5 = V'l - V'2 



(28) 



as shown in Figure [5l Now by construction 




FIG. 5: Constructing the differential measurement 3. 

l^out = 2'-^out 

and (in view of the evolution of the null geodesies about the deflector) 

lA,;„ = — 2Ai„ 

SO that 

(5 = 2 lAin- 

For an emitter at Ue = we have the rigorous result 

S = 2(/)(0) - TT. 



(29) 
(30) 
(31) 
(32) 



Of particular note is the fact that for this configuration since Aout cancels out, A plays no role, and further, the usual 
approximate argument gives the rigorously correct answer up to numerical evaluation. 

More generally, when the associated u^ are not equal, but we continue to consider opposed sources and set Ue = 0, 

5 = (01 (0) + MO) - tt) + ixiiuj - X2{uJ) (33) 

where the 0(0) are given by p^ and the x by ([^0]) . Whereas A does not enter the first term, it does not cancel from 
the second except in the limit iMe = 2Us which reproduces p2p. 

When the associated u^ are not equal, and the sources are not opposed, the differential deflection becomes 

6 = (MO) - MO)) + (Xi(^ie) " M^m)) (34) 

where again the 0(0) are given by (|12p and the x by (PO)) and we have taken Ue = 0. Whereas A does not enter the 
first term, it does not cancel from the second except in the limit lUs = 2 Ms which gives 6 = (we have one source). 



VI. DISCUSSION 

The standard argument for explaining the measured value of the deflection of light is to set the emitter and observer 
at M = so that the total deflection for a single source for grazing incidence with the Sun is given by (20(0) — 7r)|„ , 
which, when approximated by first order corrections to the null geodesic equation, gives 4u ("). For wide-angle 
deflections a similar argument leads to ([27|) . On the basis of the mathematically exact solution to the classical 
bending problem 1 have shown that these are excellent approximations. However, as regards explaining the measured 
deflection of light, these approximate arguments fall short in the sense that the bending prior to and after an apse 









TABLE III 


: «j, variable, 


n^u^ . 






rE/i?o 


xO~ 


- ^/2 (") 


x(") 


Aout (") 


A„, (") 


i'i") 


*(") 


1 


0.27 


-958.744 


959.619 


0.875246 


0.875255 


1.75050 


1.75049 


1.5 


0.40 


-1438.85 


1439.43 


0.583489 


0.583503 


1.16699 


1.16699 


1.75 


0.47 


-1678.84 


1679.34 


0.500128 


0.500145 


1.00027 


1.00027 


2 


0.53 


-1918.82 


1919.26 


0.437608 


0.437627 


0.875235 


0.875233 


2.5 


0.67 


-2398.74 


2399.09 


0.350078 


0.350101 


0.700179 


0.700178 


3 


0.80 


-2878.65 


2878.94 


0.291723 


0.291751 


0.583474 


0.583474 


3.5 


0.93 


-3358.55 


3358.80 


0.250039 


0.250072 


0.500111 


0.500111 


4 


1.07 


-3838.46 


3838.68 


0.218775 


0.218813 


0.437588 


0.437588 


4.5 


1.20 


-4318.39 


4318.58 


0.194458 


0.194501 


0.388959 


0.388958 


5 


1.33 


-4798.33 


4798.50 


0.175003 


0.175050 


0.350054 


0.350053 


5.5 


1.47 


-5278.29 


5278.45 


0.159085 


0.159137 


0.318222 


0.318221 


6 


1.60 


-5758.29 


5758.43 


0.145819 


0.145875 


0.291694 


0.291694 


6.5 


1.73 


-6238.31 


6238.44 


0.134593 


0.134593 


0.269247 


0.269247 


7 


1.87 


-6718.36 


6718.49 


0.124970 


0.125036 


0.250006 


0.250006 


8 


2.13 


-7678.58 


7678.69 


0.109331 


0.109407 


0.218737 


0.218737 


9 


2.40 


-8638.95 


8639.05 


0.0971650 


0.0972503 


0.194415 


0.194415 


10 


2.67 


-9599.52 


9599.60 


0.0874304 


0.0875252 


0.174956 


0.174956 


15 


4.00 


-14405.9 


14405.9 


0.0582079 


0.0583501 


0.116558 


0.116558 


20 


5.34 


-19220.0 


19220.1 


0.0435727 


0.0437626 


0.0873353 


0.0873353 


30 


8.02 


-28882.7 


28882.7 


0.0288895 


0.0291751 


0.0580646 


0.0580646 


40 


10.7 


-38609.6 


38609.6 


0.0214991 


0.0218813 


0.0433804 


0.0433804 


50 


13.5 


-48424.3 


48424.3 


0.0170248 


0.0175050 


0.0345299 


0.0345299 


100 


27.7 


-99811.3 


99811.3 


0.00774762 


0.00875252 


0.0165001 


0.0165001 


125 


35.6 


-128013 


128013 


0.00569625 


0.00700201 


0.0126983 


0.0126983 


150 


44.3 


-159317 


159317 


0.00417929 


0.00583501 


0.0100143 


0.0100143 


200 


68.5 


-246628 


246628 


0.00160335 


0.00437626 


0.00597961 


0.00597961 


210 


77.7 


-279667 


279667 


0.000888921 


0.00416787 


0.00505679 


0.00505679 


214.95 


90 


-323867 


323867 


0.000002635 


0.00407197 


0.00407460 


0.00407460 



are fundamentally different. The post apsidal bending involves the process of measurement and is dominated by the 
cancelation of two terms, one of which is not considered in the usual approximate argument. Whereas these more 
involved arguments may not be numerically important, conceptually they are. 

Aside form these details, the present discussion differs from others in the sense that analytic approximations to the 
exact solution ^ is not central to the discussion, as is the usual case, but rather ^ is viewed as an elementary function 
in the sense that it can be trivially approximated numerically to arbitrary accuracy using modern computational 
platforms. However, approximations have been included, relegated to the Appendices, as these serve to amplify and 
"explain" the results obtained. 
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APPENDIX A: APPROXIMATIONS TO Ai 

If 

Ue < Wj, 



then from ^ it follows that 



A,„~</)(0)-|. 



If Tj, ^ m we can take a Taylor series of (/)(0) — 7r/2 about u^ — Q which gives 

,„ „ /^ ^ 15 \ 2 /^61 15 \ o / ,^ 3465 \ , /7783 3465 \ . ,,,, 

21397 310695 \ . / 544045 765765 \ , / 400353 530675145 \ o 



12 512 y "^ V 56 256 y "^ V 8 32768 

/ 1094345069 350975625 \ „ / 3274477761 61238992815 \ ,„ .^, ^^. 
[—^2 iog^'^j"^ +( 224^+ 131072 " j "- +^("^)' 

and it is a trivial matter to go to higher orders. The first six terms in (jAip have been given previously 20]. In many 
applications the first term alone represents an adequate approximation. 

APPENDIX B: APPROXIMATIONS TO % 

Prom © and (HOI) it follows that 

cosMuyf = %2_-,,3_, (Bl) 

where 

A.^ (B2) 

and u corresponds to the observer. Equation (jBip is exact. Let us suppose that r^ 3> 2m so that 

u^ - 2u^ ~ u^ , (B3) 

and further, let us suppose that for the observer r ^ 2m, so that 

-u^ + 2u^~-u^. (B4) 



Then (IBll) takes the form 



If A = we arrive at the obvious result 



If A > let us assume that 



2 2 

11 — ?/ 

cos(xN)2 ^ ^^^. (B5) 

m1 — A 



sin(x(zi)) =^ — . (B6) 



9m'^A < 1, (B7) 

that is, 27A ^ 1. Then, for the cosmological horizon Ti, we obtain Q 

u^ ~ A (B8) 



10 

so that we can rewrite (jB5l) in the form 

? ? 

U — 11 

cos(x(u))2 ~ -^ -. (B9) 

Under the assumption r^ <^ r-j-i, that is 

ul > u^, (BIO) 

we again obtain the now somewhat less obvious result (jB6|l . 
To gauge the effect of A consider 

which reduces to 

Xa.o-Xa^o- (B12) 

Putting this into a crude cosmological context, let us take the following fiducial values: A ~ a 10~^^ cm~^, m ~ 
I (3 fO^s cm, (3 in units of lO^^ ]\^.j^^ Wj; ~ ^ 10"^ 7 in units of 100 Kpc, and m ~ ^ 10"^, S in units of Gpc. Here 

a,/3, 7 and 6 are taken to be of order unity. Now (jB3p requires - <C 10^, and (|B4|1 requires j <C 10^. For (jB7p we 
require 2a/32 <C 10^^ and for (JBIOP 807^ <C 27 10^. Relation (|B9p should therefore be a reasonable approximation in 
this context. Note that we can not set u ^ u-h as this would require 8a6^ ^ 270. Gathering this all together, we 
have the following approximation. 



Xa.o - X.^0 ^ I ^&^\ + 0{}?) ^ J — + O(A^), (B13) 

where the last term assumes j ^ 10^. For unit values of a, 7 and b then Xa=o ~ Xa^o ~ .3 which, in principle, is a 
measurable difference. However, since we can not set u ^ u-h, the cosmological context has to be refined [2l| . 

APPENDIX C: APPROXIMATIONS TO Ao„t 

If we start with the assumption r^, <C r, where r corresponds to the observer, then 

Uj, > M (CI) 

and we again obtain (/)(u) ~ (^(0) and also x — so that 

V^ ~ 20(0) - TT, (C2) 

as one finds in most texts. However, the approximation (|C1[) is inadequate. Rather, with 

u^>u (C3) 

and Uj, <'^ 10"'* to the 6 figure accuracy reported in this paper we find 



'^(«) - I + X{u) ~ l^ul - u2. (C4) 

This "explains" the almost perfect cancelation of the two terms on the left and the remarkable accuracy of the Shapiro 
- Ward formula (P7|) . It is only in the limit u = that one recovers (IC2p . 
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